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Summary 

Qh I A problem of an interface crack between two semi-planes made out of different 

. materials under an action of an in-plane loading of general tensile-shear type is treated 

in a semi-analytical manner with the help of Dirichlet-to-Neumann mappings. The 
boundaries of the crack and the interface between semi-planes are subjected to a 
curvature-dependent surface tension. The resulting system of six singular integro- 
differential equations is reduced to the system of three Fredholm equations. It is shown 
that the introduction of the curvature-dependent surface tension eliminates both 
classical integrable power singularity of the order 1/2 and an oscillating singularity 
present in a classical linear elasticity solutions. The numerical results are obtained by 
solving the original system of singular integro-differential equations by approximating 
unknown functions with Taylor polynomials. 



(N 

m 

£ — ■ 1. Introduction 



The problem of an interface crack is classical for linear elasticity and has generated an 
enormous amount of literature. The first results in this area have been obtained by Williams 
(1) on the example of a semi- infinite crack in dissimilar media. It has been shown for the 
first time that in addition to the classical singularity of the order 1/2 the stresses possess 
an oscillatory singularity of the type sin(elnr) or cos(elnr), where r is the radial distance 
from the tip of the crack. This singularity leads to the physically impossible wrinkling and 
interpenetration of the banks of the crack near the crack tips. This abnormal behavior has 
been further confirmed in the studies by England (2), Erdogan (3), Malyshev and Salganik 
(4), Rice and Sih (5) and many others. It has also been shown that the zones of overlapping 
are relatively small (10~ 7 — 1CP 4 of the crack length) in the tensile field (3), however, for 
the shear loading these zones may be comparable with the size of the crack (6). 

Different approaches have been proposed to overcome the unphysical interpenetration 
of the crack edges. One of the most well-known models has been proposed by Comninou 
(7)-(ll). The main idea is in the introduction of the contact zones near the crack tips. The 
stress singularity in this case is of the square root type r -1 / 2 and the oscillating singularity 
disappears. Interface crack problems with contact zones have been also considered in (12)- 
(14) and others. Later, Sinclair (15) introduced a crack opening angle in the framework of 
the Comninou model. 

Another approach has been proposed by Atkinson (16). In this model the interface 

Q. Jl Mech. Appl. Math. () (), 1-20 © Oxford University Press 



2 



A. Y. Zemlyanova 



is replaced by a thin strip of finite thickness. The crack is either placed inside of the 
homogeneous strip or is on the boundary of the interface layer and the parameters of the 
interface layer change continuously. Thus, in both cases the oscillating singularities do not 
occur. 

Assuming that the interfaces are rough, Mak et al (17) introduced a no-slip zone near the 
crack tips. Their solution does not contain oscillating singularities at the crack tips as well. 
Unlike the fracture model considered by Comninou where a relative shear slip is possible, 
in this model the relative slip is prohibited and the introduction of the Comninou contact 
zones is not necessary. In the paper by Boniface and Simha (18) the attempt to suppress 
an oscillating singularity has been made by choosing a specific crack opening angle which 
depends on the parameters of interacting dissimilar materials. In another recent study (19) 
the authors eliminated oscillatory singularity by introducing a surface elasticity on the crack 
boundaries. 

In this paper, we extend the curvature-dependent surface tension approach, first 
introduced in (20), (21) on the example of a straight non-interface crack, to the case of 
an interface crack between two isotropic semi-planes made out of different materials under 
an in-plane loading of general type. It is assumed that the equations of linear elasticity 
are valid everywhere in the bulk in the semi-planes. The classical linear elasticity condition 
on the crack boundaries and on the interface between semi-planes is augmented with a 
curvature-dependent surface tension. This problem in a similar formulation has been first 
approached in (22). However, this study was incomplete since an essential investigation 
of the behavior of the stresses and displacements at the crack tips and numerical results 
were absent. These results are provided in the current paper. Also it is necessary to note, 
that the regularization of the system of the singular intcgro-differential equations has been 
treated in a completely different manner from (22). The results presented in the current 
paper show that the introduction of the surface tension on the boundary of an interface 
crack leads to a complete elimination of both the singularity of the order 1/2 and the 
oscillating singularity typical for interface cracks. Thus, the introduction of surface tension 
removes classical paradoxes of linear elasticity. It is shown that one of the stresses and 
one of the derivatives of the displacements at the crack tips remains completely bounded, 
while other parameters assume logarithmic and square-logarithmic singularities. Numerical 
results presented here are achieved by approximation of the unknown functions by Taylor 
polynomials. These results are further verified by spline approximations. 

2. Statement of the problem 

Consider an upper and a lower semi-planes made out of different materials (fig. 1). The 
semi-planes are joined perfectly along x-axis with the exception of one area of debonding 
(interface crack) occupying the segment |a;| < I along the interface. The materials of the 
upper and the lower semi-planes have shear moduli /ii, /12 and Poisson's ratios v\, v-i 
correspondingly. It is assumed that the constitutive behavior of the material of each of the 
semi-planes in the bulk can be modeled with the standard equations of linear elasticity and 
only the boundary conditions are changed to incorporate the curvature dependent surface 
tension. 

Given in-plane stresses cr^, cr^, t™ x and a^, a^, r~ 2 are applied at infinity of each of 
the semi-planes. The rotations at infinity are and uj%° correspondingly. It is well-known 
that the stresses at infinity in both planes can not be independent of each other and are 
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Fig. 1 An interface crack between two dissimilar materials. 



connected by the following formulas (5): 
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Here and further the lower index "1" denotes the parameters belonging to the upper semi- 
plane and the lower index "2" denotes correspondingly the parameters of the lower semi- 
plane. At the upper and the lower boundaries of the crack the tensile f ± (x) and the shear 
g ± (x) stresses are applied. It is assumed that these stresses are in equilibrium: 

J {f + (x)+ig+(x))dx- j (f-(x)+ig-(x))dx = 0. 

Consider the jump momentum balance condition on the boundaries of the semi-planes 
(21): 

grad (c) 7 + 27ffn+[T]n = 0. (2.1) 

Here 7 is a surface tension on the boundary of the material, T denotes the Cauchy stress 
tensor, n is the unit normal to the fracture surface £ pointing into the bulk of the material, 
H = — idiv(£)n is the mean curvature, grad^-j denotes the surface gradient, and the double 
brackets [. . .] denote the jump of the quantity enclosed across the boundary of the semi- 
plane. 

Following the approach proposed in (21) for a straight non-interface crack, assume that 
the curvature dependent surface tension acts on the crack boundaries and on the interface 
between two different materials and has the following form: 



7 ± = 7^+7^div (c) n, \x\<l, 
f -7o + 7jdiv (c) n, \x\ > I, 



(2.2) 
(2.3) 



where 7o , 7i , 7oj 7i are re& l constants, the upper indexes "+", "— " and "i" describe the 
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parameters on the upper bank of the crack, on the lower bank of the crack and on the 
interface between the materials outside of the crack correspondingly. 

Assuming that the derivatives of the displacements ui,i, w 2j i and their derivatives are 
small, compute the surface divergence div^) of the inner normal n (23): 

div( C) n = u 2 ,u(x). (2.4) 

Substituting the conditions (2.2), (2.3) into the condition (2.1) and taking into account 
(2.4), one obtains the following boundary conditions: 

a+(x,0)=rfu+ lu (x,0) + f+(x), \x\<l, (2.5) 

022 fa, 0) = -Jo ut,u(x,0) + g + (x), \x\ < I, 
on the upper boundary of the crack; 

^(Z' ) = -li u 2,m( x ^) + f~( x ), \x\<l, (2.6) 

^22(^0) = lo u 2Ai( x ^) +9~{x), \x\ < I, 
on the lower boundary of the crack; and 

012(^0) - 0^0^°) = 7iW2,m (a;, 0), |a;|>Z, (2.7) 
022(^,0) - ^(^.O) = -7o«2,ii(a;,0), |a;| > I, 

on the interface between two semi-planes. Here and further "+" and "— " signs describe the 
limiting values of the stresses and the displacements on the boundaries of the semi-planes 
from the side of the upper and the lower semi-planes. 

Observe that due to the perfect attachment of the semi-planes along the interface the 
displacements m, 112, and hence their derivatives, are equal on the interface. Thus, these 
functions can be taken without "±" signs: 

uf ^x, 0) = Ui ^x, 0) = 1*1,1(2;, 0), uJi( x > 0) = u^iix, 0) = U2,i(x, 0), |a;| > I. 
3. Reduction to the system of integro-difFerential equations 

The Dirichlet-to-Neumann transforms will be used to solve the problem. It has been shown 
(21) that the Dirichlet-to-Neumann mappings for the upper semi-plane have a form: 

ai 2 (z,O) = aiU2,i(a;,O) + 0i f°° Ul ' lMdr ; (3.1) 

J -00 r-x 

CT22 (x,0) = -ai U i, ( x,0) + A r^S±. 

J -co r-x 

In a similar way as in (21) it is possible to derive that the Dirichlet-to-Neumann mappings 
for the lower semi-plane are given by the formulas: 

ai 2 (x,0)= a2U2 .i(x,0)-/3 2 ^ Ml ' lM)dr ; (3.2) 

J-oo r-x 
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(r, 0)dr 



where ct\, a 2 , (3\ and (5 2 are given constants which depend on the material parameters of 
the semi-planes: 

2fij a2 = 2f4 /?i = 2/^1+2^) = 2 M2 (A 2 +2/z 2 ) 
Ai+3/ii' A 2 + 3^ 2 ' (Ai+3^i)7r ' (A 2 + 3^ 2 )7r 

and "1 = 2(17^77' "2 = 2(A7fe7- 

Since the formulas (3.1), (3.2) are derived under the assumption that the far field loading 
in the semi-planes vanishes at infinity, it is necessary to separate the homogeneous distant 
loading held generated by the stresses applied at infinity. This leads to the following 
modifications to the boundary conditions (2.5), (2.6) and (2.7): 

a°+(x, 0) + r = 7i + ^n(^, 0) + f+(x), \x\ < I, (3.3) 

<t°+(x 7 0) + <7 = -f+u° 2 + 11 (x,0)+g + {x), \x\ < I, 

<j° 12 -(x,0)+t = -7r<m(z,0) + /-(z), \x\ < I, (3.4) 

<r 22 (x,0) + cr = To^ai^ ) + 9~( x ), M < l - 

a°+(x,0)-a° 12 (x,0)=YiU 2>ni (x,0), \x\ > I, (3.5) 

4+(x,Q)-<j° 22 -(x,0) = -7S4nM), \x\ > I, 

where all the parameters with the upper index "0" denote the stresses and the displacements 
in the absence of the loading at infinity. 
Introduce new unknown functions: 

Vi(ar) =u° 1 + 1 {x,0), ? 2 {x) = u^x, 0), \x\ < I, <p(x) = u^x, 0), \x\>l, (3.6) 

iln{x) =u° 2 + 1 {x,0), ip 2 {x) =u° 2 ; 1 {x,0), \x\<l, tl>(x) = v% A (x,0), \x\ > I. 

Substituting the Dirichlet-to-Neumann mappings (3.1), (3.2) into the boundary 
conditions (3.3), (3.4) and (3.5) leads to the following system of six singular integro- 
diffcrcntial equations for six unknown functions (3.6): 



aWi(z)+/?i/ f^L+T = 1 t^{x)+f+(x) 1 \x\<l, 

J\r\<i r-x J M>1 r-x 

-amW+hf ^^ + ( 7 = -^^ 1 (x)+ g +(x), \x\<l, 

ilrkl r-x J M>1 r-x 



i|r|</ r - a; Ji r i 



r — a; 



-^2{x)-h\ t^L^pJ i&L + a = lMx)+g -( x) , N<Z , (3.7) 

•/|r|<J r-O; J|T-|>i r-z 
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(ai - a 2 )ip(x) + 0i 

J\r 



tpi(r)dr 



:/ 

J\r\ 



ip 2 (r)dr 



-(ai - a 2 )ip(x) + /3i 



|r|</ r-a; ^|r|</ r-a; 

= -y{ip"(x), \x\ > I, 

ipi(r)dr f ip 2 (r)dr 
r P2 / 

\r\<l r-X J\ r y 



I 

J\r\>l 



+ 02 



ip(r)dr 
r — x 



ip(r)dr 



(01 + 02) [ 

J\r\>i r-x 



i\ r \<i r-x 
= -loi>'{x), \x\ > I. 

Observe that the unknown functions ipi(x), f 2 (x) and <p(x) can be easily eliminated from 
the system of the integral equations (3.7). To achieve this integrate the second, the fourth 
and the sixth equations of the system (3.7) with the Cauchy kernel and change the order of 
integration using the Poincare-Bertrand formula (24): 



J\r 



ip\(r)dr 



+Pi I In 

J\r\>l 



\ r \<i r-x 

(l-r)(l + x) 



Pm 2 i>i(x) + 0! 
ip(r)dr 



In 



r\<l 



(l-r){l + x) 



I 

J\r 

:/ 

J\r 



(l-x)(l + r) 
+ 

ip 2 (r)dr 



r — x 



+ 7o + 



(l-x)(l + r) 
tp[ (r)dr 



ipi (r)dr 



a In 



r — x 
I — x 



(3.8) 



g + (r)dr 

r\<l r — X 



\r\<l r - x 
x\ < I, 



I + X 



-a 2 

l\r\<l r-X J\r\<l 



(l-r)(l + x) 



-02 / In 

'|r|>« 



(l-r)(l + x) 



(l-x){l + r) 
+ 



tp(r)dr 



L 



r — x 
9~{r)dr 
r\<l r-x 



7o 



(l-x)(l + r) 
ip' 2 (r)dr _ 



4> 2 (r)dr 



In 



r — x 
I — x 



(3.9) 



-(«!-<*)/ ^ 

J\r\>l r- 

-02 [ 

J\r\< 



)dr 



(0! + 2 )n 2 iP(x) - [ In 

x J\r\<l 



\r\<l r-x 

x\ < I, 

(l-r)(l + x) 



I + x 



{l-x)(l + r) 



ipi(r)dr 



r — x 



(3.10) 



(l-r){l + x) 



ip 2 (r)dr 



(l-x)(l + r) 
+7o 



r — x 



- (01 + 02 



) / ^ 

J\r\>l 



(l-r)(l + x) 



(l-x)(l + r) 



ip(r)dr 



r — x 



j 

J\r\ 



ip'(r)dr 



= 0, Id > I. 



i\r\<l r-x 

Observe that for |a;| > I the formulas (3.8) and (3.9) look exactly the same with the exception 
that the terms —(3iir 2 ipi(x) and 2 n 2 ip 2 (x) drop out. Similarly, for \x\ < I the formula (3.10) 
is lacking the term — + f3 2 )ir 2 ip(x). These are the consequences of the fact that the order 
of integrals can be interchanged if at least one of the kernels is regular. 

Substituting the formulas (3.8), (3.9) and (3.10) into the first, the third and the fifth 
equations of the system (3.7) leads to the following system of three integro-differcntial 
equations: 



-7M'W+"i i- 



0W 



1p!( X ) + 



lt01 



tp'i(r)dr YoPi 



"i J|r|<; r-x a!-a 2 J M>l 



tl)'(r)dr 
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ai(ai - a 2 ) J\ r \ <t 



/ 

J\r\ 



hi 



(l-r)(l + x) 



V>i(r)<ir /3i/3 : 



r — x 



(Z-a;)(Z + r) 

A(A«2 + /?2«i) 
ai(ai-a 2 ) ,/| r |>; 

g + (r)<ir 
|r|<; r-a; 



ai 



- "2 J|r|<i 



(Z-r)(Z + a;) 



(Z-a;)(i + r) 



ip2(r)dr 



r — x 



hi 



(l-r)(l + x) 



(l-x)(l + r) 



ip(r)dr 



r — x 



Pi 



+ f + (x)-a— In 



Z — a; 



T, |x| < Z, 



7l -^) + « 2 1- 



7 #2 



Z + x 
ip' 2 (r)dr 



Oil — OL2 J\ r \<l 



In 



(l-r)(l + x) 



(l-x)(l + r) 



&2 J\ r \<l r-X 

ipi(r)dr 02 a i 



CX-l - «2 J\ r \>l 



ip'(r)dr 



r — x 



a 2 (ai — a.2) J\ r \<i 



In 



(i-r)(i + a;) 



(/-a;)(i + r) 



ip2(r)dr 



r — x 



+ 



/32(/3ic*2 +/3 2 ai) 



a 2 (ai-a 2 ) / r | >; 

#2 /" 9~(r)dr 
"2 A 



In 



(i-r)(/ + ar) 



(Z-a;)(Z + r) 
Z — x 



\p(r)dr 



r — x 



(3.11) 



Z + a; 



r, a; < 



-lW\x) + (a, - a 2 ) 1 



7o 02 
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ip'(r)dr 
r\>l r — x 
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ipi{r)dr 
r — x 

ip2 (r)dr 
r — x 

tp(r)dr 
r — x 

I — x 

In 



Z + x 



\x\ > I, 



l\r\<l r - x "2 J\ r \<i r-x 

assuming that a\ 7^ a 2 . 

Observe that the system (3.11) is not a system of singular integro-differential equations 
but rather the system of integro-differential equations with logarithmic singularities. 

Finally, to completely regularize the system in the spirit of (25), introduce the following 
new unknown functions: 



Xi(x) = ip"( x ), Xz(x) = ip2( x )i X(x) = ip"{x). 
By integrating the formulas (3.12), one obtains: 



(3.12) 



[{x) = uj (r,x)xi{r)dr + d, 

J\r\<l 
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where 



ipi(x) = / u)i(r,x)xi(r)dr + C\x + C 2 , 

J\r\<l 

ip' 2 {x) = oj (r 1 x)x2(r)dr + C 3 , 

J\r\<l 

tp 2 (x) = / ui{r 1 x)x2{r)dr + C 3 x + C 4 , 
J|r|<; 

V-'(x) = / x (r)dr, x < -I, (3.13) 

r 

ip(x) = / (x — r)x(r)dr, x < —I, 

J —oo 

ip'(x) = - x(r)dr, x > I, 

J X 

ip(x) = — (x — r)x(r)dr, x > I, 

J X 

w {r,x) = |^ ^ \-lx] vi(r,x) = (x-r)w (r,x). 



The representations (3.13) contain four real constants Ci, C2, C3 and C4. These 
constants will be found from additional physical conditions which will be stated later. The 
corresponding integration constants for the function x( x ) are chosen to be equal to zero to 
guarantee the appropriate decay of the function tjj(x) at infinity and the existence of all the 
integrals in (3.11). 

Finally, substituting the formulas (3.12), (3.13) into the system (3.11) reduces this system 
to the final system of Fredholm equations the theory of which is well-developed. 



4. Additional conditions 

The representations (3.13) contain four real constants C\, C 2 , C3 and C4. To fix these 
constants four additional real conditions coming from the physics of the problem will be 
stated. Firstly, it is reasonable to assume that the total stresses applied to the boundaries 
of the crack must be equal to zero, which gives two real conditions: 

J ((rt 2 -(^ 22 )dx = 0, (4.1) 

/ <Tu)dx = 0. (4.2) 

The equation (4.1) can be expressed through the unknown functions and the constants Ca- 
using the boundary conditions (3.3), (3.4): 

J {l^'M +%i>' 2 {x) -g + {x) +g-(x))dx = 0. 
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At the same time to rewrite the condition (4.2) in terms of the unknown functions and the 
constants Cj it is better to use the Dirichlet-to-Neumann mappings (3.1), (3.2). 

To make the solution physically reasonable it is necessary for the displacements to be 
single-valued along the boundary of the crack. In other words, if the position of the left 
end of the crack is fixed then the displacement of the right end of the crack must be the 
same along both the upper and the lower semi-planes. This leads to the following two real 
conditions: 

y «i -ui A )dx = 0, (4.3) 

/ Ki -«2,i)<fc = 0. (4.4) 

These conditions are easy to rewrite in terms of the unknown functions and the constants 
Cj as : 

J (ipi(x) - (p 2 (x))dx = 0, 

J (M*)- Mx))dx = 0. 

The conditions (4.1)-(4.4) allow us to fix all the constants in the formulas (3.13). 



5. Singularities of the stresses and derivatives of the displacements at the 
crack tips 

Consider the system of singular intcgro-differential equations (3.9). Observe that the 
functions <pi(x), y>i{x), tp(x), ipi(x), ip2{x) and tp(x) can not have integrable power 
singularities of any power (including pure imaginary powers) at the crack tips x = ±1. 
Indeed, assume that the singularities are possible and consider temporarily only the point 
x = I. Then 

<p k (x) = A k (x-iy + o((x-iy), (5.1) 

i>'l{x) = B k {x-iy + ((x-iy), 

where the index k is equal to 1, 2 or missing and A k , B k are real coefficients. It follows 
then that the functions ip' k {x) 7 "4> k (x) are bounded at the point x — I. From the second, 
the fourth and the sixth equations of the system (3.9) it follows then that the functions 
<p k (x) may have at most logarithmic singularities. Using this fact and the first, the third 
and the fifth equations of the system (3.9), deduce then that the functions il> k (x) may have 
at most logarithmic or square-logarithmic singularities. Thus, we have a contradiction with 
the initial representations (5.1). 

Hence, the functions Pk(x), "4> k (x) can not have integrable power singularities of any 
power. This means, that unlike in the linear theory of elastic fracture mechanics (LEFM), 
the stresses and the derivatives of the displacements do not possess singularities of the order 
1/2 or oscillating singularities with a pure imaginary power. Thus, the introduction of the 
surface tension on the boundary of the crack removes these important contradictions of the 
linear theory. 
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On the other hand, it can be seen from the system (3.9) that the functions <fik(x), 4> k (x) 
can have logarithmic singularities of the following form: 

Vk(x)=A k \n{x-l) + 0{l), (5.2) 

V4'(x) = B k ln 2 (x -l) + C k ln(x - I) + 0(1), 

where A k , B k and C k are real coefficients. Then, as before it follows that the functions 
ip' k (x), ip k (x) are bounded at the point x = I. Hence, using the boundary conditions (3.3), 
(3.4), (3.5) it can be seen that the stresses are bounded at the point x = I and the 
stresses af 2 have a logarithmic singularity of the same type as the functions ip k (x) (5.2). 
The derivative of the displacement uf 1 is bounded and the derivative uf 1 has the singularity 
of the same type as <p k (x) . Obviously, the same conclusions are valid at the left tip of the 
crack x = —I. Thus, the singularities of the stresses and of the displacements are of the 
same type as for a curvilinear non-interface crack (26) and hence, the surface tension model 
does not differentiate between an interface and a non-interface cracks. Observe also, that 
for the conclusion of this section to stand it is essential that the coefficients 7^ and are 
all non-zero. 



6. Numerical results 

To produce the numerical results for the stresses and the displacements, the singular 
integro-differential system (3.7) which has a significantly more convenient form than the 
regularized system (3.11) will be solved numerically. Two different approaches will be 
considered. The first approach consists in representing the unknown functions by Taylor 
polynomials, substituting these polynomials into the system (3.7) and consequently solving 
the resulting system of linear algebraic equations. This approach has been found to be the 
most convenient and to produce the best numerical results. For a verification purposes only 
this approach has been corroborated by representing the unknown functions with splines. 

First of all, the unknown functions (p k (x) , i\) k (x) are approximated by Taylor polynomials 
with unknown coefficients: 

N N N+l 

= ^ a \x k , ^(i)=^a^, ^(x)=^afx- fe , (6.1) 

fc=0 fc=0 k=l 

N+2 N+2 N+l 



X 



Mx) = E b l xk > Mx) = E b t xk > ^) = E h l 

k=0 k=0 fe=l 

The idea behind the representations (6.1) is that all of the series have the same number of 
terms (N + 1) except for the series for ipi(x) and tp2(x) which have two additional terms 
which correspond to the constants Ci, C2, C3 and C4. In addition, the functions <p{x), 
ip(x) must decline to zero at infinity which guarantees the convergence of all the integrals 
involved. 

Substituting the representations (6.1) into the system (3.7) leads to the following integrals 
which are easy to compute analytically: 

L ^ L = _ y ( IJ ^ 

i„i w r — x L — 1 ] — k 
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/ ?*-± M<l, (6.2) 

J\r\>ir-x J + k 

[ — —E 1 "^*-^ M>'> 

|r|>i r * j=ljj¥fc 3 k 

Assume that the tensile and the shear stresses applied to the boundaries of the crack have 
the following Taylor series representations: 

/ + (z) = E/ fe + * fe , n*) = E ( 6 - 3 ) 

fc=0 fc=0 
N N 

k=0 k=0 

By substituting the formulas (6.1), (6.2), (6.3) into the system (3.7) and making sure 
that the coefficients by the powers of x are the same on both sides of the equations, obtain 
the following system of linear algebraic equations: 

- 7l +(fc + 2)(fc + l)b\ +2 = -TS ok + /+, k = 0, . . . , N, 
+ 7o +(fc + = + ,g+, fc = 0, . . . , N, 

Q!2& fc -^2 - 2^ T— a 3 l + E -r T V 

+ 7 f(fc + 2)(fc + l)&2 +2 = _ T( 5 0fe + k = 0,...,N, 

-a 2 at - A ( - E 4^^^ + E ^^"^ I ^ 
-7(7 ( fc + l ) b l+i = -vSok+gk, k = 0,...,N, 
(a 1 -a 2 )6l + (/3l+ ^ ) E ^^^-AE 1 !^^ 

3=1,3^ ' J=0 
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'fa E 1 -,y +fc ffl fc+J ' - 7 i (fc - 2)(fc - 1)63fe - 2 = °' k = l,...,N + l, 



1 _ (_i\j+k 

'fa E ? + fy J+k fo(k l)6ti = 0, fc = 1, . . . , iV + 1, 

where (5^ is the Kronecker symbol. 

It is also possible to rewrite the additional conditions (4.1)-(4.4) in terms of the Taylor 
series representations (6.1). Thus, the following four additional equations for the coefficients 
a*, bj must be satisfied: 

N+2 N+2 

7+ J2 h \e k {l (-l) fe ) + 7o" E b l e "^ (-!)") = °' ( 6 - 5 ) 

fe=0 fe=0 
N+2 N+2 

7l + E b\ke k -\l {-If- 1 ) + 7 f E tike'-Hl M)^ 1 ) = 0, (6-6) 
fe=i fe=i 

E°^ i (1 " fc ( ;T +1) -E^ fc+i (1 "i;T +1) = °> ^ 

fc=0 fc=0 

E^ii^-E^^^P-. ,,S, 

fc=0 k=0 

Hence, the equations (6.4)-(6.8) form a system of QN + 10 real linear algebraic 
equations with 6N + 10 real unknowns. Solving this system gives the Taylor polynomial 
approximations (6.1) of the unknown functions. 

The results obtained by this method have been verified on several examples by 
additionally solving the system (3.7) numerically using spline approximations of the 
unknown functions. To obtain these approximations, first of all, rewrite the system (3.7) in 
the form which involves integration over finite intervals only: 



<Pi{r)dr f ip (r)dr 

r P 1 / rl 

\ r \<i r-x J H<1 rx-V 



aiMx) + ft + ft + T = lM{x) + ;+(,), W < lt 



J\r\<l r-x J M< i rx-V 

a*h{x)-(h[ ^L-ftl ^ + r = - 7l -^) + /"(*), N<i, 
J\r\<i r-x J\r\<i rx - I 



-a 2 <p 2 (x)-fa[ ^^-/3 2 / i^T)± + a = ^' 2 (x)+g-(x), \x\ < I, (6.9) 
J\r\<i r-x J M<1 rx-V 
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{ai - a 2 )M x ) + ^ I ^ + ^ I ^ + (ft + ft) / 

J\r\<i rx-P J ]rl<l rx-P J M 



(Pojr)dr 

l\r\<r'X-l- J\ r \<l ix - i- j| r | <; r-X 

= 1 \l- 4 (2x 2 i, {x) + 4x 3 iP' (x) + xVo'W), \x\ < I, 
, \ i \ n i2 f ipi(r)dr . j2 f ip2(r)dr f ip (r)dr 

J\r\<l rX - 1 J\r\<l rx - L J\r\<l r ~ x 

= f r 2 {xM x ) + x2 ^{x)), \x\ < I, 

where ip (x) = £<p ) , M x ) = §V> (£) ■ 

Take the linear spline approximations of the functions ifi(x), y>i(x) and fo{x): 

TV XV 

tpk{x) = -j^(x-Xj-i) + ^—{xj-x), x e[xj-i,Xj], j = l,...,2N, fc = 0,l,2, (6.10) 

where the points = — I + jl/N, j = 0, . . . , 2N, generate the partition of the segment 
[—1, I], and the values of the functions (fik{xj) at these points are equal to w k . To guarantee 
the proper decline of the function <p(x) at infinity (and, hence, the existence of all the 
corresponding integrals) it is necessary that Wq — 0. Thus, the approximations (6.10) 
contain 6iV + 2 unknowns w k . 

To be able to differentiate the approximations of the functions ipi(x), tp2(x) and ipo(x) 
twice it is necessary to take the higher order approximations of these functions which can 
be achieved by approximating the functions tpi(x), ip2(x) and ipo(x) by cubic splines (27): 

M x )=S{{x), ie[iM,4 j = l,...2N, (6.11) 
S J k (x) = x^f + ^-( Xj x) 3 (6.12) 

^zij (x - xj-!) + - fa - x )- 

It is easy to see that the functions ipk{x) represented by (6.11) and their second derivatives 
ip'k(x) are continuous through the points Xj, j = 1, . . . , (2N — 1). To achieve the continuity 
of the derivatives ip' k (x) at these points we must satisfy the following conditions 

z{- 1 + Az{ + 4 +1 = ^(yi- 1 2y{ + y{ +1 ), j = 0,...,2N. (6.13) 

Also assume that in the formulas (6.12), (6.13) 

^ = *r +1 = o, y~ k l = vl N+l = o, fc = o,i,2. 

Similarly, to guarantee an adequate decline of the function tjj(x) at infinity it is necessary 
to set Dq = 0. Thus, the formulas (6.11), (6.12) contain 127V + 5 unknowns y 3 k , z k and the 
equations (6.13) provide 6N + 3 conditions to find these unknowns. 

Substituting the representations (6.10), (6.11) into the system (6.9) and into the 
additional conditions (4.1)-(4.4) produces the remaining 12iV + 4 linear equations. Observe, 
that the equations of the system (6.9) are satisfied numerically only at the 2N points 
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Fig. 2 Comparison of the results obtained by the Taylor polynomials and the spline 
approximations. 



£j = {xj_i + Xj)/2, j = 1,...,2N, which are the midpoints of the segments [xj-i,Xj]. 
Solving the resulting system of linear algebraic equations numerically results in the spline 
approximations of the functions ip^x), ipk( x ) at the points Xj, j = 0, . . . , 2N. 

Comparison of the results obtained by the approximations with the Taylor polynomials 
and with splines is presented on the fig. 2. The results have been obtained for \i\ = = 70, 
vi=i/2 = 0.3, I = 1, and the following surface tension constants: 7^ = 7;f = 0.01, 7^ = 
7(7 = 0.01, 7] = 7g = 0.005. The construction is subjected to the pure tensile loading a = 1 
at infinity and the boundaries of the crack are free from external stresses. The solid line 
depicts the shape of the crack obtained with Taylor polynomial approximations for N = 30, 
the dotted line represents the results obtained with Taylor polynomial approximations for 
N = 50, and, finally, the diamond marks represent the spline approximations for N = 30. 
It is easy to see that all the graphs are almost indistinguishable from each other. Since the 
Taylor polynomial approximations were, in author's opinion, more computationally efficient, 
the results presented further in this paper were obtained by this method. 

The results for the interface crack problem in the case of the LEFM are well-known and 
have been derived by many authors for different types of configurations and loading. In 
this paper, we compare the results obtained with the curvature-dependent surface tension 
model with the classical results of England (2). In the case of the absent loading at infinity 
and a uniform pressure f + {x) = f~(x) = 0, g + (x) = g~(x) = T, applied to the boundaries 
of the crack, the displacements of the boundaries of the crack are derived in an explicit form 
according to the following formulas (2): 




(6.14) 
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Fig. 3 Influence of the surface tension on the shape of the crack. 



u 2 (x) 



2^ 2 (l + a) 



(I 2 - x 2 ) 1/2 cos 7 In 



+ x 



\x\ < I, 



(6.15) 



where 



Mi + M2«i 1 , 

a= , 7= — In a. 

M2+Ml K l 27T 



l+x 



ill 



Observe the presence of the oscillating singularities given by the terms cos ^7 In 

the formulas (6.14), (6.15). 

The influence of the magnitude of the surface tension on the shape of the crack is presented 
on the fig. 3. The shape of the crack is plotted for I = 1, /ii = 70, /J2 = 80, v\ = 0.3, 
V2 = 0.35, and the following surface tension constants: 7i" = 7f = 7i = 7tf = 7t7 = 7o = 7- 
The constant 7 takes on the following five values: 7 = (corresponds to the classical LEFM 
solution obtained by the formulas (6.14), (6.15)), 7 = 0.001, 7 = 0.01, 7 = 0.1 and 7 = 1.0. 
Observe, that the introduction of the surface tension principally changes the behavior of 
the stresses and the displacements at the crack tips, however, the magnitude of the surface 
tension constants has little influence on the displacements of the boundaries of the crack. 
Observe also, that the displacements are almost equal for 7 = 0.001, 7 = 0.01 and 7 = 0.1, 
and only for 7 = 1.0 a noticeable difference can be seen. 

The shape of the crack and the stresses on the crack boundary are shown on the fig. 4 
for different types of loading of the construction. The computations are made for = 70, 
A*2 = 80, vi = 0.3, v 2 = 0.35, I = 1, 7^ = 0.02, 7f = 0.01, 7^ = 0.01, % = 0.02, 7J = 0.01, 
7q = —0.01. The graphs (a), (b) correspond to the shape of the crack and the stresses 
on the crack boundary for the case when the construction is subjected to the pure tensile 
loading a = 1 at infinity in the absence of any other loading; the graphs (c), (d) to the case 
of the pure shear loading at infinity t = 1; the graphs (c), (f) to the case of the loading 
,f + (x) = f~(x) = 0, g + (x) = g~{x) = \ll 2 — x 2 applied on the boundary of the crack and 
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Fig. 4 Shape of the crack and stresses on the crack boundary for different types of loading. 

the absence of the loading at infinity. Note that the edges of the crack on the graph (c) 
self- intersect. This is because of the appearance of the contact zones on the boundaries of 
the crack in the case of the shear loading. The self-intersection may be removed if these 
zones are explicitly accounted for during the solution of the problem. Also observe that 
with the exception of the stresses af 2 on the fig. 4(f), the stresses on the boundaries of 
the crack are close to zero everywhere except for the small area near the tips of the crack. 
This is due to the fact that the curvature dependent surface tension is small in magnitude 
everywhere except for the small zones near the crack tips. However, presence of this tension 
essentially changes the character of the singularities at the crack tips. 

The dependence of the maximal stresses and the singularity coefficients in the construction 
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Fig. 5 Dependence of the maximal stresses and the singularity coefficients on the Poisson's ratio 
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Fig. 6 Dependence of the maximal stresses and the singularity coefficients on the shear modulus 
Hi. 
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Fig. 7 Dependence of the maximal stresses and the singularity coefficients on the parameter 7. 
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of the Poisson's ratio of one of the semi-planes v\ is shown of the fig. 5. The other parameters 
of the construction are \i\ = [i 2 = 70, v 2 = 0.30, I — 1, 7^ = 7f = 7^" = 7^ = 0.01, 
71 = 7o = 0.005. The figure (a) corresponds to the pure tensile loading at infinity u = 1 
and the figure (b) corresponds to the pure shear loading r = 1. The coefficients fci and k 2 
are the coefficients by the square logarithmic singularities of the stresses a 12 : 

a+ 2 {x) = kt -x) + OQn(l - x)), a^(x) = k 2 In (I - x) + 0(ln(Z - x)). 

Observe that the maximums of the stresses are achieved at the tips of the crack. Similar 
graphs for the dependence of the maximal stresses and the singularity coefficients on the 
shear modulus \i\ are shown on the fig. 6. All the parameters are the same as before with 
the exception of the changing parameter /zi, and v\ — 0.3. The graphs of the dependence of 
the maximal stresses and the singularity coefficients on the parameter 7, 7^ = 7f = 7^ = 
7o~ = 7, 7i = 7o = 0.005 are shown on the fig. 7. The other parameters of the construction 
are \x\ — 70, ji 2 = 80, v\ = 0.3, v 2 = 0.35, / = 1. The figure (a) as before correspond to the 
pure tensile loading a = 1 at infinity and the figure (b) to the pure shear loading r = 1. 

Conclusions 

In this paper, a new model of brittle fracture first proposed in (23) and later used to study 
fracture problems in (20)-(22) is applied to the problem of an interface fracture between 
two dissimilar isotropic semi-planes under the tensile-shear in-plane loading. The classical 
boundary condition of linear elasticity is augmented with a curvature-dependent surface 
tension. 

The physical problem is reduced to the system of singular integro-differential equations 
by applying Dirichlet-to-Neumann mappings. The system of singular integro-differential 
equations is further reduced to the system of Fredholm integral equations. It is shown that 
unlike in the classical solutions of linear elasticity, the introduction of the surface tension 
leads to the solutions which do not possess integrable power singularities of the order 1/2 
or the oscillating singularities, thus making the solutions more practically feasible. The 
integrable singularities of the logarithmic type may still be present. The most important 
problem arising from the presented solution is in developing more general models of surface 
tension which may allow to eliminate logarithmic singularities as well, thus making the 
solutions completely bounded in the neighborhoods of the crack tips. This is the subject of 
ongoing research. 

The numerical solution of the system is obtained by approximating the unknown functions 
by Taylor polynomials and is further verified by spline collocation method. The results of 
computations are presented for different types of mechanical parameters. 
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